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We develop a theory of nonlinear cosmological perturbations on superhorizon scales for generic 
single-field inflation. Our inflaton is described by the Lagrangian of the form W(X, cj>) — G(X, </>)□</> 
with X = —d tl (l>di_ l (j)/2, which is no longer equivalent to a perfect fluid. This model is more general 
than k-inflation, and is called G-inflation. A general nonlinear solution for the metric and the 
scalar field is obtained at second order in gradient expansion. We derive a simple master equation 
governing the large-scale evolution of the nonlinear curvature perturbation. It turns out that the 
nonlinear evolution equation is deduced as a straightforward extension of the corresponding linear 
equation for the curvature perturbation on uniform <j> hypersurfaces. 

PACS numbers: 98.80.-k, 98.90.Cq 



I. INTRODUCTION 



Non-Gaussianity in primordial fluctuations is one of the most powerful tools to distinguish different models of 
inflation (see, e.g., Ref. [l[ and references therein). To quantify the amount of non-Gaussianity and clarify non- 
' Gaussian observational signatures, it is important to develop methods to deal with nonlinear cosmological pertur- 
O" 1 bations. Second-order perturbation theory 043 is frequently used for this purpose. Spatial gradient expansion is 
employed as well in the literature [oT-|27j. The former can be used to describe the generation and evolution of pri- 
bJ)' mordial perturbations inside the horizon, while the latter can deal with the classical nonlinear evolution after horizon 
exit. Therefore, it is important to develop both methods and to use them complementarily. 

In this paper, we are interested in the classical evolution of nonlinear cosmological perturbations on superhorizon 
scales. This has been addressed extensively in the context of the separate universe approach [l3| or, equivalent ly, 
the SN formalism [1, [M CLl G1 ■ The <5A^ formalism is the zeroth-order truncation of gradient expansion. However, 
t"***" ' higher-order contributions in gradient expansion can be important for extracting more detailed information about 
CO . non-Gaussianity from primordial fluctuations. Indeed, it has been argued that, in the presence of a large effect from 
the decaying mode due to slow-roll violation, the second-order corrections play a crucial role in the superhorizon 
t— I , evolution of the curvature perturbation already at linear order (28rl3fj| . as well as at nonlinear order [H, [23|. In order 
' to study time evolution of the curvature perturbation on superhorizon scales in the context of non-Gaussianity it 
CO | i s necessary to develop nonlinear theory of cosmological perturbations valid up to second order in spatial gradient 
I ■ expansion. 

The gradient expansion technique has been applied up to second order to a universe dominated by a canonical scalar 
field and by a perfect fluid with a constant equation of state parameter, P/p =const 0, |2(J. The formulae have 
been extended to be capable of a universe filled with a non-canonical scalar field described by a generic Lagrangian 
of the form W(—d fJ, 4>d fJl 4>/2,(j>), as well as a universe dominated by a perfect fluid with a general equation of state, 
P = P(p) [HI, Those systems are characterized by a single scalar degree of freedom, and hence one expects that 
a single master variable governs the evolution of scalar perturbations even at nonlinear order. By virtue of gradient 
expansion, one can indeed derive a simple evolution equation for an appropriately defined master variable £K^ L : 

^L" + 2 -^ L ' + ^ (2) i?[^ L ] =G(e 4 ), (1.1) 
z 4 

where the prime represents differentiation with respect to the conformal time, e is the small expansion parameter, 
and the other quantities will be defined in the rest of the paper. This equation is to be compared with its linear 
counterpart: 

n^" + 2-Kl i »'-c 2 s A1ll i » = 0, (1.2) 

z 
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from which one notices the correspondence between the linear and nonlinear evolution equations. Gradient expansion 
can be applied also to a multi-component system, yielding the formalism "beyond SN" developed in a recent paper [25j . 

The purpose of this paper is to extend the gradient expansion formalism further to include a more general class 
of scalar-field theories obeying a second-order equation of motion. The scalar-field Lagrangian we consider is of the 
form W(— d fJ, (/>d l j l (j)/2, </>) — G(—d ll (f)d l j r (f)/2,(j))D(f). Inflation driven by this scalar field is more general than k-inflation 
and is called G-inflation. It is known that k-inflation [HI, HH with the Lagrangian W{— d^^d^j 2, <fi) is equivalently 
described by a perfect fluid cosmology. However, in the presence of GO(f>, the scalar field is no longer equivalent to a 
perfect fluid and behaves as a imperfect fluid [33l |3~H ] . The authors of Ref. (26j investigated superhorizon conservation 
of the curvature perturbation from G-inflation at zeroth order in gradient expansion, and Gao worked out the zcroth- 
order analysis in the context of the most general single-field inflation model [27j . In this paper, we present a general 
analytic solution for the metric and the scalar field for G-inflation at second order in gradient expansion. By doing so 
we extend the previous result for a perfect fluid (2~H and show that the nonlinear evolution equation of the form 
is deduced straightforwardly from the corresponding linear result even in the case of G-inflation. 

This paper is organized as follows. In the next section, we define the non-canonical scalar-field theory that we 
consider in this paper. In Sec. IIII| we develop a theory of nonlinear cosmological perturbations on superhorizon scales 
and derive the field equations employing gradient expansion. We then integrate the relevant field equations to obtain 
a general solution for the metric and the scalar field in Sec. IIV1 The issue of defining appropriately the nonlinear 
curvature perturbation is addressed and the evolution equation for that variable is derived in Sec. [V] Section I VII is 
devoted to a summary of this paper and discussion. 



II. G-INFLATION 



In this paper we study a generic inflation model driven by a single scalar field. We go beyond k-inflation for which 
the Lagrangian for the scalar field is given by an arbitrary function of 4> and X := —g^ u d^(j)d v (t>/2, C = W(X, (j>), and 
consider the scalar field described by 

I = J d i xy/^g[W(X,(f>)-G(X,(t)p(f ) ], (2.1) 

where G is also an arbitrary function of 4> and X. We assume that <fi is minimally coupled to gravity. Although the 
above action depends upon the second derivative of <fi through □</> = (^V^V^, the resulting field equations for <f> 
and g^u remain second order. In this sense the above action gives rise to a more general single-field inflation model 
than k-inflation, i.e., G-inflation |35l-l37l|. The same scalar- field Lagrangian is used in the context of dark energy and 
called kinet ic g ravity braiding [331] . Interesting cosmological applications of the Lagrangian (|2.ip can also be found, 
e.g., in j38l443| . In fact, the most general inflation model with second-order field equations was proposed in [37j 
based on Horndcski's scalar-tensor theory [4^, EHJ . However, in this paper we focus on the action (|2.1j) which belongs 
to a subclass of the most general single-field inflation model, because it involves sufficiently new and interesting 
ingredients while avoiding unwanted complexity. Throughout the paper we use Planck units, M p \ = 1, and assume 
that the vector —g^V ' v <\> is timelike and future-directed. (The assumption is reasonable because we are not interested 
in a too inhomogeneous universe.) 

The equation of motion for (f> is given by 

V M [(W x -G$- G x n0) W - GxVl] + W,/, - G^U<j> = 0, (2.2) 

where the subscripts X and <f> stand for differentiation with respect to X and (/>, respectively. More explicitly, we have 

WxUcf) - Wxx{^,NuW<t>V v <t>) - 2W^ X X + W <t> - 2(G - G^ X X)U$ 

+G X [(V„V^)(VV^ - (n^) 2 + ii M „ WV>] + 2G^(V M V„</>)(VW0) + 2G 00 A 

-G X x{VV x <j> - ^ A n0)(V M V»V„</>V A </> = 0. (2.3) 

The energy-momentum tensor of the scalar field is given by 

T„ u = WxVpipVut + Wg^ - (V M GV„0 + V„GV M 0) + <? M „V A GV A </) - Gx^V^V^. (2.4) 

It is well known that k-inflation allows for an equivalent description in terms of a perfect fluid, i.e., the energy- 
momentum tensor reduces to that of a perfect fluid with a four-velocity = d tl (f>/\/2X . However, as emphasized 
in [33| , for Gx 7^ the energy-momentum tensor cannot be expressed in a perfect-fluid form in general. This imperfect 
nature characterizes the crucial difference between G- and k-inflation. 
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III. NONLINEAR COSMOLOGICAL PERTURBATIONS 



In this section we shall develop a theory of nonlinear cosmological perturbations on superhorizon scales, following 
Ref. plj . For this purpose we employ the Arnowitt-Deser-Misner (ADM) formalism and perform a gradient expansion 
in the uniform expansion slicing and the time-slicc-orthogonal threading. 

A. The ADM decomposition 

Employing the (3 + l)-decomposition of the metric, we write 

ds 2 = g^dx^dx" = -a 2 dt 2 + %(dx i + fi l dt)(dx 3 + (5 3 dt), (3.1) 

where a is the lapse function and (5 l is the shift vector. Here, latin indices run over 1,2,3. We introduce the unit 
vector normal to the constant t hypersurfaces, 

n dx" = -adt, n^d, =-(d t - 0*$). (3.2) 

a 

The extrinsic curvature of constant t hypersurfaces is given by 

Ku = ViUj = (d t % - DiPj - . ( 3 - 3 ) 

where Di is the covariant derivative associated with the spatial metric 7^. The spatial metric and the extrinsic 
curvature can further be expressed in a convenient form as 

% = o a (t)e a *^7tf, (3-4) 
Ka = a 2 (t)e 2 ^ (±K 7ij + A t ^j , (3.5) 

where a(t) is the scale factor of a fiducial homogeneous Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime, 
the determinant of 7^ is normalised to unity, det 7^ = 1, and A^ is the trace- free part of Kij, 7' y Aij = 0. The trace 
K := 7 lJ is explicitly written as 



K=i 

a 



3{H + d t V) - Di? , (3.6) 



where H = H(t) is the Hubble parameter defined by H := d\na(t)/dt. In deriving Eq. p.6[) dt(det%j) = 7 l ° dtfij = 
was used. Hereafter, in order to simplify the equations we choose the spatial coordinates appropriately to set 

/3 ? = 0. (3.7) 

We call this choice of spatial coordinates as the time-slice-orthogonal threading. 

With fa = all the independent components of the energy-momentum tensor (|2.4p are expressed as 



E 
-Ji 



T^n v = (W x - G x V4>){d ± <t>y -W- d ± Gd ± 4> - fidiGdjcf), (3.8) 
l^n" = [{Wx - G x u<i>)di<t> - d.G] - d±G8i<f>, (3.9) 
Ty, (3.10) 



where d± :— n^d^. 



Let us now move on to the (3 + l)-decomposition of the Einstein equations. In the ADM Language, the Einstein 
equations are separated into four constraints (the Hamiltonian constraint and three momentum constraints) and six 
dynamical equations for the spatial metric. The constraints are 



_1 o„, , 2 

a 
2 



2 i?[e 2 ^ 7 ] + -K 2 - A Z3 A l i = 2E, (3.11) 



8 t K - e-^Dj (e 3 ^) = J,, (3.12) 
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where i?[e 2 ^7] is the Ricci scalar constructed from the metric e 2 ^"fij and Di is the covariant derivative with respect 
to 7^ . The spatial indices here are raised or lowered by and 7^, respectively. As for the dynamical equations, 
the following first-order equations for the spatial metric (ip, 7^) are deduced from the definition of the extrinsic 
curvature p.3[) : 

dxi> = + (3-13) 
a 3 

d^ij = 2A ij . (3.14) 
The dynamical equations for the extrinsic curvature (K, A^j) are given by 



K 2 D 2 a 1 

d±K = -— - WAh + — - - (E + 3P) , (3.15) 



1 



TF 1 r„ , o„, , „ iTF 



d±Aij = -KA ij +2A i k A kj + - D^a - [ify [e 2 ^ 7 ] - S y ] , (3.16) 



P := ^-"e-^Sij, (3.17) 



o 

where 

1 

3 

D 2 := 7 y DiDj, and i?^ [e 2 ^] is the Ricci tensor constructed from the metric e 2 ^jij. The trace-free projection 
operator [. . .] is defined for an arbitrary tensor as 

[Q l3 ] TF := Qa - pyyUQu- (3-18) 

For the purpose of solving the Einstein equations, the most convenient choice of the temporal coordinate is such 
that the expansion K is uniform and takes the form: 

K(t,x) = 3H(t). (3.19) 

Hereafter we call this gauge choice with (|3.7[) the uniform expansion gauge. Adopting this gauge choice, Eq. (j3.13[) 
reduces simply to 

dt4) = H(a-l)=:H6a(t,x). (3.20) 

From this, if we take the uniform expansion gauge, we can see that the time evolution of the curvature perturbation ip is 
caused by the inhomogeneous part of the lapse function Sa(t, x l ) only. It is related to the non-adiabatic perturbation. 



B. Gradient expansion: basic assumptions and the order of various terms 



In the gradient expansion approach, we introduce a flat FLRW universe characterized by (a(t), <f>o(t)) as a back- 
ground and suppose that the characteristic length scale L = a/k, where k is a wavenumber of a perturbation, is longer 
than the Hubble length scale 1/H of the background, HL ^> 1. We use e := 1/(HL) = k/(aH) as a small parameter 
to keep track of the order of various terms and expand the equations in terms of e, so that a spatial derivative acting 
on a perturbation raises the order by e. 

The background flat FLRW universe characterized by (a(t), <j>o(t)) satisfies the Einstein equations, 



H 2 (t) 



1 



;Po, 



H(t) = -±(p + P ), 



and the scalar-field equation of motion, 



Jo + 3H Jo = (W4, - 2XG H - 2G^ x Xq> 



(3.21) 



(3.22) 



Here, an overdot (') denotes differentiation with respect to t, and a subscript indicates the corresponding background 
quantity, i.e., Wq := W(X , <f>o), (Wx)o ■— Wx(Xo, <fo), etc., where X := </> 2 /2. The background energy density and 
pressure, po and Pq, are given by 



Po - 
Po 



-W + 2X(W X 
W-2X(Gx4>- 



- 3HG X ^ 
Ga 1 



Ga 



(3.23) 
(3.24) 
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while J7o is defined as 



Jo 



(3.25) 



Wx^-ZG^ + QHGxX . 
L Jo 

If VF and G do not depend on </>, the right hand side of Eq. (|3.22[) vanishes and hence J$ is conserved. In this case, J7b 

is the Noether current associated with the shift symmetry <p — > (f> + c. Note that the above quantities may be written 
in a different way as 

Po = Jofo ~W + 2(XG^) , (3.26) 

Po + Po = Jofo ~ 2(G X X0) O . (3.27) 
Note also that the scalar-field equation of motion (|3.22p can be written as 

G4> + 39 Jo + £4 = 0, (3.28) 

where 

g(t):=Sx-3Q(Gxi>)o, (3-29) 

0(t) := H - (G x Xj>) , (3.30) 

£*(t) := [2W X ^ - + GHG^xXj) - 2XG U ^ ^ , (3.31) 
£ x (t) : = \w x + 2X1V XX + 9i?G x + &HG X xi>o - 2G^ - 2XG <px ] q ■ (3.32) 

These functions will also be used later. 

Since the background FLRW universe must be recovered at zeroth order in gradient expansion, the spatial metric 
must take the locally homogeneous and isotropic form in the limit e — > 0. This leads to the following assumption: 

dt-ra = C(e 2 ). (3.33) 

This assumption is justified as follows [3, Il9l - [25| . If dtjij were 0(e), the leading term would correspond to ho- 
mogeneous and anisotropic perturbations, which are known to decay quickly. We may therefore reasonably assume 
dtjij = 0(e 2 ) and not dtjij = 0(e). However, ip and 7^ (without any derivatives acting on them) are of order 0(1). 

Using the assumption (|3.33[) made above and the basic equations derived in the previous subsection, one can now 
deduce the order of various terms in gradient expansion. First, from Eq. Q3.14p we see that 

A y =0(e 2 ). (3.34) 

Substituting Eq. (|3.34[) into Eq. (|3.12[) under the gauge condition (|3.19[) . we obtain Ji — 0(e 3 ). Then, this condition 
combined with the definition (|3.9[) implies that di8cf> = 0(e 3 ), where 5<p(t,x) := </>(i,x) — 4>o(t). The same equations 
also imply that diG = 0(e 3 ). By absorbing a homogeneous part of d(j> into 4>o (and redefining a(t) accordingly), we 
have 

6<p = (D(e 2 ). (3.35) 
It is clear from the condition (13.341) and the Hamiltonian constraint p. lip that 

SE :=E{t,x l ) -p (t) =0(e 2 ). (3.36) 
Since the definition Q3.8P tells that E — p = max{0(<5</>), 0(d±(f> — <9 t 0o)}, we see d t (S(f>) — tftoSa = 0(e 2 ), leading to 

Sa = 0(e 2 ). (3.37) 

Then, it follows immediately from Eq. p.20p that 

d t ip = 0(e 2 ). (3.38) 
Similarly, for the spatial energy-momentum component we find 

5P := P(t, x) - P (t) = 0(e 2 ). (3.39) 
In summary, we have evaluated the order of various quantities as follows: 
^ = 0(1), 7«=0(1), 

5a = 0(e 2 ), 5<j) = 0(e 2 ), 6E = 0(e 2 ), 6P = 0(e 2 ), Ay = 0(e 2 ), 

d aij = 0(e 2 ), d t ^ = 0{e 2 ), /3 8 =0(e 3 ), d t G = 0(e 3 ), [5 y ] TF = 0(e 6 ), (3.40) 
where the assumptions made have also been included. 



C. Field equations up to 0(e 2 ) in gradient expansion 

Keeping the order of various terms (|3.40[) in mind, let us derive the governing equations in the uniform expansion 
gauge. The Hamiltonian and momentum constraints are 



i?[e 2V, 7 ] = 2SE + (D(e i ), 

The evolution equations for the spatial metric are given by 

d t i> = HSa + 0(e 4 ), d tll] = 2A l} + 0(e 4 ), 
while the evolution equations for the extrinsic curvature are 

dtAii = -ZHA l0 - [i? y -[e^ 7 ]] TF + 0(e 4 ), 

-d t H = -3H 2 -l(E + 3P)+0(e 4 ). 
a 2 

Note that with the help of the background equations Eq. p.45[) can be recast into 

The components of the energy-momentum tensor are expanded as 

E = 2XW X — W + 6HG x dj_<j) - 2XG + C(e 6 ), 
P= W -d ± Gd ± <p + 0{e 6 ), 
—Ji = Jodi{6ct>) - (G x <j>) di(8X) + 0(e 5 ), 

where 

X = (d ± cf > ) 2 /2 + 0(e e ), 

SX := X - X = 4> Q d t (54>) - 2X 5a + 0(e 4 ). 

Finally, noting that D<f> = —d\(j> — 3Hd±(f> + 0(e 4 ), the scalar-field equation of motion (|2.3[) reduces to 

W x {dl4> + 3Hd ± <j>) + 2XW xx dl(b + 2W 4>X X -W+- 2(G - G^ x X){dl<t> + 3Hd ± 4>) 
+6G x [d ± (HX) + 3H 2 X] - AXC^d 2 ^ - 2G^X + 6HG xx Xd ± X = C(e 4 ). 

This equation can also be written in a slightly simpler form as 

d±J + 3HJ = W tf> - 2XG^ - 2XG 4>x d 2 ^ + C(e 4 ), 

where 

J = W x d ± <p - 2G 4> d ± (t) + QHG X X. 

It can be seen that Eq. p.53[) takes exactly the same form as the background scalar-field equation of motion 
under the identification dt o d±. Now Eq. (|3.47p can be written using J as 

E = Jd x <f>-W + 2XG^. 

From Eqs. (pH)3")) and (f3~551) we find 

d ± E = -3H(E + P)+ C(e 4 ). 

This equation is nothing but the conservation law, n M V V T^ V = 0. 
Combining Eq. (|3.45[) with Eq. (|3.56[) . we obtain 

d t [a 2 (5E)} = 0(e 4 ). 



3.41) 
3.42) 

3.43) 

3.44) 
3.45) 

3.46) 

3.47) 
3.48) 
3.49) 

3.50) 
3.51) 

3.52) 
3.53) 
3.54) 



3.55) 
3.56) 

3.57) 
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Wc can expand Eq. (|3.47|) in terms of 6<j) and 8X, and thus 8E can be expressed as 

SE = £^{t)8<j> + £ x (t)SX + 0(e 4 ). 



(3.58) 



This equation relates 8<fi and 8X with a solution to the simple equation p.57|) . With the help of Eq. (|3.51[) . Eq. (|3.58j) 
can be regarded as an equation relating 5(f) and 8a. Similarly, one can express 5P as 



SP = 



Using Eq. (|3.46[) . one has 



^d t {a 3 [j {S(l>) - (Gx4)o(SX)] } - ( Po + P )(Sa) + 0(e 4 ). 



3 

d t {a 3 [Jo(54>) - (G X <P) (8X)\ } = -^-SE + 0(e 4 ), 



(3.59) 



(3.60) 



which can easily be integrated once to give another independent equation relating 8<p and 8a. In the next section, we 
will give a general solution to the above set of equations. 



IV. GENERAL SOLUTION 



Having thus derived all the relevant equations up to second order in gradient expansion, let us now present a general 
solution. First, since ip = 0(1) and d t ip = 0(e 2 ), we find 



i!j= (0) C' / '(x)+O(e 2 ), 



(4.1) 



where ^C^(x) is an integration constant which is an arbitrary function of the spatial coordinates x. Here and 
hereafter, the superscript (n) indicates that the quantity is of order e™. Similarly, it follows from 7^ = 0(1) and 
d t7ij =0{e 2 ) that 



7u = (°)C^x) + 0( e 2 ), 



(4.2) 



where (°)C^(x) is a 3 x 3 matrix with a unit determinant whose components depend only on the spatial coordinates. 
The evolution equations (|3.43j) can then be integrated to determine the 0(e 2 ) terms in tp and 7^ as 



(4.3) 
(4.4) 



where t* is some initial time and integration constants of 0(e 2 ) have been absorbed to (°'C^(x) and ^ - , C^ (x). 
Now Eq. (|3.44|) can be integrated to give 



,0 = (0)c*(x) + J H(t r )Sa(t', X )dt' + 0(e 4 ), 
7« = (O) 05(x) + 2 f Aijif »df + 0(e 4 ), 



A, 



a 3 (t) 



( 2 )^(x) y a(i')dt' + (2) ^(x) 



0(e 4 ), 



(4.5) 



where 



1 



(2)^.. _ 1(2)^(0)^^ + (d^C^d^C* - ( )A (0) £>/ 0) C^) 



(4.6) 



Here, ( o )0T fc * is the inverse matrix of (°)07., (2) %(x) := fly [(°) C"?] and ( 2 >i?(x) := fl[(°)C"T] arc the Ricci tensor 
and the Ricci scalar constructed from the zeroth-order spatial metric ^°^C^ (x), and (°>D is the covariant derivative 
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associated with '°'C^. Note that WCT-?' ( 2 ).F]y = by dehnition. The integration constant, ( 2 )C^(x), is a sym- 
metric matrix whose components depend only on the spatial coordinates and which satisfies the traceless condition 
(o) cl ij (2) C A = q Substituting the above result to Eq. JO}, we arrive at 



* df 



t. « 3 (*') 

Next, it is straightforward to integrate Eq. f|3 . 5T[) to obtain 



a(t")dt" + (2) C4(x) 



df 



0(e 4 



SE 



1 



« 2 W 



.(2) 



/C(x)+0(e 4 ), 



(4.7) 



(4.8) 



where ( 2 )/C(x) is an arbitrary function of the spatial coordinates. With this solution for SE, Eqs. (|3.58[) and (|3.60[) 
reduce to 



£^t)dd> + £ x (t)SX = -L-( 2 )/C(x) + 0(e 4 ), 
a 2 (t) 



0, 



{a 3 [MS® - (G x i)o(6X)]} 



1 



.(2) 



and the latter equation can further be integrated to give 

( 2 >C*(x) 



JoW>) - (Gx4>)o(SX) 



3a 2 (t) 

(2)/C(x) 



/C(x)+0(e 4 ), 



a 3 (i) 3a 3 (t) 



a(i')^' + ^(e 4 ), 



(4.9) 
(4.10) 

(4.11) 



where we have introduced another integration constant ^C x (x). 

With the help of Eqs. (|3.5ip . one can solve the system of equations (|4.9[) and (|4.11[) , leading to 



1 

A 



(G 



and 



S a = d t - 



1 

2X S(t) 



3a(t) 



(2) £ , gxft) (2) . 



a 2 (t) a 3 (t) 



+0(e 4 ), 



( 2 >/C 39(t) (a , 
a 2 (i) a 3 (i) 



0(e 4 



(4.12) 



(4.13) 



where .4 := {E^Gx'P + £x»7)lo- I n deriving the above solution we have used the background scalar-field equation of 
motion (|3~28| . 

Finally, substituting Eq. (|4.13l) to Eq. (|4.3j) . we obtain 



^ = (0)c*+ / dt'H(t')d, 



t» 



2X g(t f ) 



e(f) 



a{t")dt'' 



a 2 (i') a 3 (t') 



(o) 



00 



2X ga 2 (f) 



2^0 
a(tf) 



a(f)M- 



( 2 )^ + ^)( 2 )c4+0( e 4 ), 



(4.14) 



where we performed integration by parts and used the background equation. 

So far we have introduced five integration constants, WC^(x), C°)Cy(x), ( 2 )Cy(x), ( 2 '/C(x), and ( 2 )C x (x), upon 
solving the field equations up to C(e 2 ). Here, it should be pointed out that they are not independent. Indeed, 
Eqs. ()3.4ip and (13.421) impose the following constraints among the integration constants: 

^x) = ^W), 

e -^O*(p )cr M0) D . \j m C*> (2) C A (x) j = Q, (2) C x (x) +0(£ 5 ); 

1 
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c^ 2 >i?(x)+0(e 5 ), 



(4.15) 
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where ^ 2 ^i?(x) := i?[e 2< >c * (°)(7 7 ] is the Ricci scalar constracted from the metric e 2< )c * ^Cy- Here, ^ 2 ^-R(x) should 
not be confused with ( 2 'i?(x). The latter is the Ricci scalar constructed from '"'C^ and not from e 2<0>c "'' i°)(j7^. 



Explicitly, 



(2) 



i?(x) = [< 2 >i?(x) - 2 (2(°)fl 2 (%* + (»)ct^,(°)c^/%*) 



(4.16) 



Note that the third equation is automatically satisfied provided that the last equation holds, as can be verified by 
using Eq. (|4~6j) . 

In summary, we have integrated the field equations up to second order in gradient expansion and obtained the 
following solution for generic single-field inflation: 



6E 



(2) 



2a 2 
1 



0(e% 



5a 



Aa 2 

(0) 



0(e 4 ), 



1 



2X Ga 2 
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a{t')dt' 



(2; 



i?(x) 36 



C^(> 



dt' 



H5<j) 
+ - 






H 


[(-; 


2X Ga 2 



2</>o 
t' 

a(t")dt' 



( 2 )C x (x) 



^M + ^( 2 )CMx) 



0(e 4 ), 



0(e 4 ), 



( 2 )^-(x) / a(i')^' + (2) ^(x) 



7« = (0) ^(x) 



(2) ^-(x) 



dt' 

t, aW) Ju 



+ 0(6 4 ), 

a(f)df + ( 2 >d(x) 



dt 1 



0(e 4 



(4.17) 



The x-dependent integration constants, (°)C" / ', ^C^, ^C x and ^C^, satisfy the following conditions: 
^. = (°'CJ„ det(C°>C£) = 1, (2 ^ = ( 2 '^, ^C™ ^ = 0, 

e -3f)C* C0) a 7ife(0 )£) . ^»C* (2) C A^ = ft (2)^ 



(4.18) 



Before closing this section, we remark that the gauge condition (|3.7p remains unchanged under a purely spatial 
coordinate transformation 



x l -> a; 4 = /*(x). (4.19) 

(°)C^ contains three residual gauge degrees of freedc 
the number of degrees of freedom associated with each integration constant is summarized as follows: 



This means that the zeroth-order spatial metric ^ contains three residual gauge degrees of freedom. Therefore, 



C'C^ •■• 1 scalar growing mode = 1 component, 

^Cij ■ ■ ■ 2 tensor growing modes = 5 components — 3 gauge, 

i 2 )(jx ... i scalar decaying mode = 1 component, 

^ Cpj ■ ■ ■ 2 tensor decaying modes = 5 components — 3 constraints. 



(4.20) 



V. NONLINEAR CURVATURE PERTURBATION 



In this section, we will define a new variable which is a nonlinear generalization of the curvature perturbation up to 
0(e 2 ) in gradient expansion. We will show that this variable satisfies a nonlinear second-order differential equation, 
and, as in Ref. f22j , the equation can be deduced as a generalization of the corresponding linear perturbation equation. 
To do so, one should notice the following fact on the definition of the curvature perturbation: in linear theory the 
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curvature perturbation is named so because it is directly related to the three-dimensional Ricci scalar; ip may be 
called so at fully nonlinear order in perturbations and at leading order in gradient expansion; and, as pointed out in 
Ref. [lij], ip is no longer appropriate to be called so at second order in gradient expansion. To define the curvature 
perturbation appropriately at 0(e 2 ), one needs to take into account the contribution from 7^-. Let us denote this 
contribution as x- We carefully define the curvature perturbation as a sum of -0 and x so that the new variable 
reproduces the correct result in the linear limit. In what follows we remove the subscript from the background 
quantities since there will be no danger of confusion. 

A. Assumptions and definitions 

As mentioned in the previous section, we still have residual spatial gauge degrees freedom, which we are going to fix 
appropriately. To do so, we assume that the contribution from gravitational waves to 7^ is negligible and consider the 
contribution from scalar-type perturbations only. We may then choose the spatial coordinates so that 7^ coincides 
with the fiat metric at sufficiently late times during inflation, 

jij — > Sij (t —> 00). (5-1) 

In reality, the limit t — > 00 may be reasonably interpreted as t — > ii ate where ii a t e is some time close to the end of 
inflation. Up to 0(e 2 ), this condition completely removes the residual three gauge degrees of freedom. 

We wish to define appropriately a nonlinear curvature perturbation on uniform </> hypersurfaces {8<j){t, x) = 0) and 
derive a nonlinear evolution equation for the perturbation. The nonlinear result on uniform <j) hypersurfaces is to 
be compared with the linear result for G-inflation [35j . However, in the previous section the general solution was 
derived in the uniform expansion gauge. For our purpose we will therefore go from the uniform expansion gauge to 
the uniform </> gauge 1 . It is clear that at leading order in gradient expansion the uniform expansion gauge coincides 
with the uniform </> gauge. In this case, it would be appropriate simply to define ip to be the nonlinear curvature 
perturbation. At second order in gradient expansion, however, this is not the correct way of defining the nonlinear 
curvature perturbation. We must extract the appropriate scalar part x from 7^ , which will yield an extra contribution 
to the total curvature perturabtion, giving a correct definition of the nonlinear curvature perturbation at 0(e 2 ). 

Let us use the subscripts K and u to indicate the quantity in the uniform K and <f> gauges, respectively, so that in 
what follows the subscript K is attached to the solution derived in the previous section. First, we derive the relation 
between ipx and ip u up to C(e 2 ). In general, one must consider a nonlinear transformation between different time 
slices. The detailed description on this issue can be found in Ref. [25|. However, thanks to the fact that 54>k = C(e 2 ), 
one can go from the uniform K gauge to the uniform </> gauge by the transformation analogous to the familiar linear 
gauge transformation. Thus, ipu is obtained as 

i> u = i>K-^S^ K +0(e 3 ). (5.2) 

One might think that the shift vector fi % u appears in this new variable as a result of the gauge transformation, but 
(3 1 can always be gauged away by using a spatial coordinate transformation. The general solution for ip u valid up to 
0(e 2 ) is thus given by the linear combination of the solution for ipK and 8cf>K displayed in Eq. ()4.17|) . Note here that 
the spatial metric 7^ remains the same at C(e 2 ) accuracy under the change from the uniform K gauge to the uniform 
<f> gauge: 

7^ = 7^ + 0(e 4 ). (5.3) 

We now turn to the issue of appropriately defining a nonlinear curvature perturbation to C(e 2 ) accuracy. Let 
us denote the linear curvature perturbation in the uniform <p gauge by 7l[f n . In the linear limit, ip reduces to the 
longitudinal component H\ ln of scalar perturbations, while x to traceless component H^ ln : 

4>^HY U \ X^H^\ (5.4) 



The gauge in which (j> is uniform is sometimes called the unitary gauge. The unitary gauge does not coincide with the comoving gauge 
in G-inflation, as emphasized in [35| . 
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The linear curvature perturbation is given by 7?. Lln = (H^ m + Hj, ln /3)Y. Here, we have followed Ref. [46[ and 
the perturbations are expanded in scalar harmonics Y satisfying [did 1 + fc 2 ) Y^ = 0, with the summation over k 
suppressed for simplicity. The spatial metric in the linear limit is expressed as 

= « 2 (% + 2Hl in Y6 l3 + 2H^ a Y ij ) , (5.5) 

where = k~ 2 [didj — (l/3)5ijdid l ] Y^. Since tp corresponds to H^ m , one can read off from the above expression 
that jij = Sij + 2HJf ul Yij in the linear limit. Thus, our task is to extract from 7^ the scalar component x that reduces 
to H^ m in linear limit. It was shown in Ref. [22| that by using the inverse Laplacian operator on the flat background, 
A -1 , one can naturally define \ as 

X ■= ^A- 1 [0* e - 3 W^( 7ij - Sii)] • (5.6) 

In terms of \ defined above, the nonlinear curvature perturbation is defined, to C(e 2 ), as 

V\ NL := V + I • (5.7) 

As is clear from Eq. (|5.6p . extracting \ generally requires a spatially nonlocal operation. However, as we will 
see in the next subsection, in the uniform <j> gauge supplemented with the asymptotic condition on the spatial 
coordinates (|5.1|) . it is possible to obtain the explicit expression for the nonlinear version of \ from our solution (|4.17[) 
without any nonlocal operation. 



B. Solution 



We start with presenting an explicit expression for ip u . It follows from Eqs. f|4. 1 T[) and (|5 .2[) that 
i' u = (°)C*(x) + < 2 )^'(x) + f R (t) ( 2 )i?(x) + f x (t) ( 2 >C*(x) + G(6 4 ). 



(5.8) 



Note here that although the integration constant ( 2 )C^(x) was absorbed into the redefinition of (°>C^(x) in the 
previous section, we do not do so in this section for later convenience. The time-dependent functions /#(£) and fc(t) 
are defined as 



/*(*) := 
fx(t) ■= 



dt' \ (P + P) 
2a 2 } 24>A 



G 



3 a 



a{t")dt" 



H 

2~xg 



e 



a{t")dt'' 



' dt' 



(p + P)£x 'ZHQ 



2<M 2X< 5. 

Since ~fij U coincides with "fij k up to 0(e 2 ), it is straightforward to see 

llJU = (°)Cg.(x) + ™C%(x) + 2 5j ,(t)< a >*i i (x) + 2 5A (t)< a >Cg(x) + 0(e 4 ), 

where 



9F(t) 



dt' 

t, W) Ju 



a(t")dt", g A (t) 



dt' 



(5.9) 
(5.10) 

(5.11) 
(5.12) 



The integration constants ^C?- and ^C?- are determined from the condition (|5.1[) as 

(0) ^ = * y , 



(5.13) 



We now have < -°- ) C^ = 5^, and hence ^Rij(x) = R^l^C] = 0. This simplifies the explicit expression for ^R(x) 

and ( 2 ).Fjj(x); they are given solely in terms of ( >C^ and the usual derivative operator c^. 
Substituting Eq. ()5.11[) to the definition (|5.6[) . we obtain 



X11 



-^-> [g F (t) - 5f(oo)] ^ M*) - ^(°°)] + °( £ )■ 



(5.14) 
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It is easy to verify that the linear limit of Xu reduces consistently to HJf ln Y. We then finally arrive at the following 
explicit solution for the appropriately defined nonlinear curvature perturbation in the uniform <f> gauge: 

^NL = (0) C ^( X ) + (2)G*(x) + < 2) i?(x) 

Let us comment on the dependence of 5H^ L on the initial fiducial time t*. One may take i* as the time when our 
nonlinear superhorizon solution is matched to the perturbative solution whose initial condition is fixed deep inside 
the horizon. Then, should not depend on the choice of i*, though apparent dependences are found in the lower 
bounds of the integrals fn,{t), f x (t), g F {t), and gA{t). Actually, in the same way as discussed in Ref. [22j|, one can 
check that 5H^ L is invariant under the infinitesimal shift — > f » + 5t„ . 



/*(*) 



g F (t) g F (oo) 



12 



12 



(2)CX(j 



fx(t) - 



9A(t) , 3a(oo) 



(5.15) 



C. Second-order differential equation 

Having obtained explicitly the solution £H^ L in Eq. (|5.15[) . now we are going to deduce the second-order differential 
equation that obeys at 0(e 2 ) accuracy. For this purpose, we rewrite /fl(t) and f x (t) in terms of 

-■ G ^ (5.16) 



(-) 



This is a generalization of familiar "z" in the Mukhanov- Sasaki equation (4?]], and reduces indeed to a^J (p + P)/Hc s 
in the case of k- inflation. With some manipulation, it is found that fn(t) and f x (t) can be rewritten as 

, If a(r/)ir/ 1 f irj ["' „ 1 f" dr/ p' 2 „ 

1 [ v dr,' p dr]' 



^—-Lim-'liwy (5 - 17) 

where the conformal time defined by dr\ = dt/a(t) was used instead of t, and rj* corresponds to the fiducial initial 
time. Further, it is convenient to express them in the form 



Mr?) = F( V *) - F(r,) - -g F (r)) , f x {rf) = -g A ( V ) + D(r,*) - D{rf), (5.18) 



where we defined 



The functions D(r]) and ^(r?) are defined so that D, F — >• as 7] —> 0. It is important to notice that -D(f?) is the 
decaying mode in the long- wavelength limit, i.e., at leading order in gradient expansion, in the linear theory, satisfying 

*' 

D" + 2—D' = 0, (5.20) 

while F{rj) is the 0(k 2 ) correction to the growing (constant) mode satisfying 

z' 

F" + 2-F' + c 2 s =0, (5.21) 

where we assume that the growing mode solution is of the form 1 + k 2 F(rj) + 0(k i ). In the above equations the 
prime stands for differentiation with respect to the conformal time and c 2 is the sound speed squared of the scalar 
fluctuations defined as 

c ' := ffp - f w : =^(- 9 ' + 0g ^)o- ( 5 - 22 ) 

Using D and F, Eq. f|5 . 1 5[) can be written as 

K L (V) = (0) C^(x) + (2 ^C m {x) - ( 2 >i?(x)F(7/) - WC x (x)D{ V ) + G(e 4 ), (5.23) 
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where time-independent terms of C(e 2 ) are collectively absorbed to ( 2 )C K (x). It turns out that the solution can be 
expressed simply in terms of the two time-dependent functions corresponding to the decaying mode and the 0(k 2 ) 
correction to the growing mode in the linear theory. This shows that, within 0(e 2 ) accuracy in gradient expansion, 
the curvature perturbation obeys the following nonlinear second-order differential equation: 

9C" + 2^ L ' + f {2) R[< L ] = ^(e 4 ) , (5.24) 

where ( 2 'i?[fH„ L ] is the Ricci scalar of the metric Sij exp (25K^ L ). Equation f|5 . 24[) is our main result. It is easy to 
see that in the linear limit Eq. (|5.24j) reproduces the previous result for the curvature perturbation in the unitary 
gauge HI], 

K^ n " + 2-7^ in ' - c 2 A7^ in = , (5.25) 

where A denotes the Laplacian operator on the flat background. 

Equation (|5.24|) can be regarded as the master equation for the nonlinear supcrhorizon curvature perturbation at 
second order in gradient expansion. It must, however, be used with caution, since it is derived under the assumption 
that the decaying mode is negligible at leading order in gradient expansion. Moreover, if one the right-hand side of 
(|5.24p set to exactly zero, this master equation becomes a closed equation, and it be a useful approximation to a full 
nonlinear solution on the Hubble horizon scales or even on scales somewhat smaller than the Hubble radius. 



VI. SUMMARY AND DISCUSSION 



In this paper, we have developed a theory of nonlinear cosmological perturbations on superhorizon scales for G- 
inflation, for which the inflaton Lagrangian is given by W(X, <j>) — G(X, </>)□</>. In the case of Gx = 0, i.e., k-inflation, 
the energy-momentum tensor for the scalar field is equivalent to that of a perfect fluid. In the case of G-inflation, 
however, it can no longer be recast into a perfect fluid form, and hence its imperfect nature shows up when the 
inhomogeneity of the Universe is considered. We have solved the field equations using spatial gradient expansion in 
terms of a small parameter e :— k/(aH), where k is a wavenumber, and obtained a general solution for the metric 
and the scalar field up to 0(e 2 ). 

We have introduced an appropriately defined variable for the nonlinear curvature perturbation in the uniform 
<f> gauge, 5H^ L . Upon linearization, this variable reduces to the previously defined linear curvature perturbation 
7i^ ln on uniform <fi hypcrsurfaccs. Then, it has been shown that £K^ L satisfies a nonlinear second-order differential 
equation (|5.24p . which is a natural extension of the linear perturbation equation for 7?-„. We believe that our result 



can further be extended to include generalized G-inflation, i.e., the most general single-field inflation model [371 ] . 
though the computation required would be much more complicated. 

The nonlinear evolution of perturbations, and hence the amount of non-Gaussianity, are affected by the C(e 2 ) 
corrections if, for example, there is a stage during which the slow-roll conditions are violated. Calculating the three 
point correlation function of curvature perturbations including the C(e 2 ) corrections will be addressed in a future 
publication. 
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